Finding Apparent Horizons in Dynamic 3D Numerical Spacetimes

We have developed a general method for finding apparent horizons in 3D numerical relativity. Instead of solving for the partial differential equation describing the location of the apparent horizons, we expand the closed 2D surfaces in terms of symmetric trace--free tensors and solve for the expansion coefficients using a minimization procedure. Our method is applied to a number of different spacetimes, including numerically constructed spacetimes containing highly distorted axisymmetric black holes in spherical coordinates, and 3D rotating, and colliding black holes in Cartesian coordinates.Comment: 19 pages, 13 figures, LaTex, to appear in Phys. Rev. D. Minor changes mad

New Coordinate Systems for Axisymmetric Black Hole Collisions

We describe a numerical grid generating procedure to construct new classes of orthogonal coordinate systems that are specially adapted to binary black hole spacetimes. The new coordinates offer an alternative approach to the conventional \v{C}ade\v{z} coordinates, in addition to providing a potentially more stable and flexible platform to extend previous calculations of binary black hole collisions.Comment: 3 pages, 5 postscript figures, LaTeX, uses mprocl.sty (available at http://shemesh.fiz.huji.ac.il/MG8/submission.html) To appear in the proceedings of the Marcel Grossmann 8 (Jerusalem, 1997

Towards the Final Fate of an Unstable Black String

Black strings, one class of higher dimensional analogues of black holes, were shown to be unstable to long wavelength perturbations by Gregory and Laflamme in 1992, via a linear analysis. We revisit the problem through numerical solution of the full equations of motion, and focus on trying to determine the end-state of a perturbed, unstable black string. Our preliminary results show that such a spacetime tends towards a solution resembling a sequence of spherical black holes connected by thin black strings, at least at intermediate times. However, our code fails then, primarily due to large gradients that develop in metric functions, as the coordinate system we use is not well adapted to the nature of the unfolding solution. We are thus unable to determine how close the solution we see is to the final end-state, though we do observe rich dynamical behavior of the system in the intermediate stages.Comment: 17 pages, 7 figure

Reflection of a Lieb-Liniger wave packet from the hard-wall potential

Nonequilibrium dynamics of a Lieb-Liniger system in the presence of the hard-wall potential is studied. We demonstrate that a time-dependent wave function, which describes quantum dynamics of a Lieb-Liniger wave packet comprised of N particles, can be found by solving an $N$-dimensional Fourier transform; this follows from the symmetry properties of the many-body eigenstates in the presence of the hard-wall potential. The presented formalism is employed to numerically calculate reflection of a few-body wave packet from the hard wall for various interaction strengths and incident momenta.Comment: revised version, improved notation, Fig. 5 adde

Tracking Black Holes in Numerical Relativity

This work addresses and solves the problem of generically tracking black hole event horizons in computational simulation of black hole interactions. Solutions of the hyperbolic eikonal equation, solved on a curved spacetime manifold containing black hole sources, are employed in development of a robust tracking method capable of continuously monitoring arbitrary changes of topology in the event horizon, as well as arbitrary numbers of gravitational sources. The method makes use of continuous families of level set viscosity solutions of the eikonal equation with identification of the black hole event horizon obtained by the signature feature of discontinuity formation in the eikonal's solution. The method is employed in the analysis of the event horizon for the asymmetric merger in a binary black hole system. In this first such three dimensional analysis, we establish both qualitative and quantitative physics for the asymmetric collision; including: 1. Bounds on the topology of the throat connecting the holes following merger, 2. Time of merger, and 3. Continuous accounting for the surface of section areas of the black hole sources.Comment: 14 pages, 16 figure

The Asymmetric Merger of Black Holes

We study event horizons of non-axisymmetric black holes and show how features found in axisymmetric studies of colliding black holes and of toroidal black holes are non-generic and how new features emerge. Most of the details of black hole formation and black hole merger are known only in the axisymmetric case, in which numerical evolution has successfully produced dynamical space-times. The work that is presented here uses a new approach to construct the geometry of the event horizon, not by locating it in a given spacetime, but by direct construction. In the axisymmetric case, our method produces the familiar pair-of-pants structure found in previous numerical simulations of black hole mergers, as well as event horizons that go through a toroidal epoch as discovered in the collapse of rotating matter. The main purpose of this paper is to show how new - substantially different - features emerge in the non-axisymmetric case. In particular, we show how black holes generically go through a toroidal phase before they become spherical, and how this fits together with the merger of black holes.Comment: 28 pages, 10 figures, uses REVTEX. Improved quality figures and additional color images are provided at http://www.phyast.pitt.edu/~shusa/EH

Adaptive Event Horizon Tracking and Critical Phenomena in Binary Black Hole Coalescence

This work establishes critical phenomena in the topological transition of black hole coalescence. We describe and validate a computational front tracking event horizon solver, developed for generic studies of the black hole coalescence problem. We then apply this to the Kastor - Traschen axisymmetric analytic solution of the extremal Maxwell - Einstein black hole merger with cosmological constant. The surprising result of this computational analysis is a power law scaling of the minimal throat proportional to time. The minimal throat connecting the two holes obeys this power law during a short time immediately at the beginning of merger. We also confirm the behavior analytically. Thus, at least in one axisymmetric situation a critical phenomenon exists. We give arguments for a broader universality class than the restricted requirements of the Kastor - Traschen solution.Comment: 13 pages, 20 figures Corrected labels on figures 17 through 20. Corrected typos in references. Added some comment

Generalized harmonic formulation in spherical symmetry

In this pedagogically structured article, we describe a generalized harmonic formulation of the Einstein equations in spherical symmetry which is regular at the origin. The generalized harmonic approach has attracted significant attention in numerical relativity over the past few years, especially as applied to the problem of binary inspiral and merger. A key issue when using the technique is the choice of the gauge source functions, and recent work has provided several prescriptions for gauge drivers designed to evolve these functions in a controlled way. We numerically investigate the parameter spaces of some of these drivers in the context of fully non-linear collapse of a real, massless scalar field, and determine nearly optimal parameter settings for specific situations. Surprisingly, we find that many of the drivers that perform well in 3+1 calculations that use Cartesian coordinates, are considerably less effective in spherical symmetry, where some of them are, in fact, unstable.Comment: 47 pages, 15 figures. v2: Minor corrections, including 2 added references; journal version

Exact Solutions for the Intrinsic Geometry of Black Hole Coalescence

We describe the null geometry of a multiple black hole event horizon in terms of a conformal rescaling of a flat space null hypersurface. For the prolate spheroidal case, we show that the method reproduces the pair-of-pants shaped horizon found in the numerical simulation of the head-on-collision of black holes. For the oblate case, it reproduces the initially toroidal event horizon found in the numerical simulation of collapse of a rotating cluster. The analytic nature of the approach makes further conclusions possible, such as a bearing on the hoop conjecture. From a time reversed point of view, the approach yields a description of the past event horizon of a fissioning white hole, which can be used as null data for the characteristic evolution of the exterior space-time.Comment: 21 pages, 6 figures, revtex, to appear in Phys. Rev.